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Abstract 



o 



Nonextensive statistics in a Blast- Wave model (TBW) is implemented to describe the identified hadron production in relativistic 
p+p and nucleus-nucleus collisions. Incorporating the core and corona components within the TBW formalism allows us to describe 
simultaneously some of the major observations in hadronic observables at the Relativistic Heavy-Ion Collider (RHIC): the Number 
of Constituent Quark Scaling (NCQ), the large radial and elliptic flow, the efifect of gluon saturation and the suppression of hadron 
^ production at high transverse momentum {pr) due to jet quenching. In this formalism, the NCQ scaling at RHIC appears as a 
l~5 consequence of non-equilibrium process. Our study also provides concise reference distributions with a least fit of the available 
[ — . experimental data for future experiments and models. 
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\ Several intriguing features were discovered in relativistic 
heavy ion collisions HI |2> S 01 when particles emerging from 
the Quark-Gluon Plasma were detected by the experiments at 
RHIC. In AuH-Au collisions, identified particle yields integrated 
pver the transverse momentum range around the center-of-mass 
rapidity window have been shown to be at equilibrium at the 
phemical freeze-out in a statistical analysis LL JJ. The hydrody- 
namic model with proper equation of state and initial condition 
can describe the anisotropic flow with small shear viscosity and 
provides the notion of "perfect liquid" [6, 7]. Furthermore, the 
^transverse momentum distributions of identified particles can 
be described in a hydrodynamic -inspired model with a compact 
set of parameters [1, 5, 8, 9, 10]. 

■ However, in the intermediate pt range, particle production 
exhibits grouping between baryons and mesons with baryons 
having relatively higher yield and larger elliptic flow than the 
mesons |Ill[Tlll. This feature of constituent quark scaling is 
not present in the hydrodynamics. A microscopic quark coa- 
lescence at the hadronization seems to be inescapable iH, [12I1 . 
At even higher pj, hard perturbative QCD processes (jets) are 
relevant. Absorption of jets in the medium formed in Ah-A 
collisions has been used for studying the properties of the 
QGP [T3I [l3l ■ Even though hydro^namics with space-time 
evolution from an initial condition [6] is so far the most real- 
istic simulation for bulk matter produced in relativistic heavy 
ion collisions, its applicabiUty is expected to breakdown for 
pH-p and peripheral A+A collisions at RHIC. Recent study 



showed that hydrodynamics can not replace the microscopic 
hadronic cascade at the late stage regardless of freeze-out and 
equation of state one chooses [15] because the particle inter- 
actions may be dominated by non-equilibrium hadronic pro- 
cesses lfl6ll . In A-hA collisions, the fluctuations at initial impact 
due to Color-Glass Condensate (CGC) formation or individual 
nucleon-nucleon collision may not be completely washed out 
by subsequent interactions in either the QGP phase or hadronic 
phase [17]. These effects leave footprints in the spectra at low 
and intermediate pj- 

With its development and success of nonextensive statistics 
(also known as Tsallis statistics [18]) in dealing with non- 
equilibrated complex systems in condensed matter, many au- 
thors have utilized Tsallis statistics to understand the particle 
production in high-energy and nuclear physics l flo[[l9[l20[ 21, 
,22, .231 ■ Although the implications and understanding of the 
consequences of such an application are still under investiga- 
tion, the usual Boltzmann distribution in an nij exponential 
form can be readily re-written as an mj power-law (Levy) func- 
tion ill: 



InrnTdrriTdy 



cc (1 + ^^^)-l/(?-l) 



(1) 
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where the left-hand side is the invariant differential particle 
yield and ^ is a parameter characterizing the degree of non- 
equilibrium, niT = -^rrP- + pj is the transverse mass of the 
given particle with mass of m and T is related to the temper- 
ature of the system. The distribution can be derived from the 
usual procedure in statistical mechanics, starting from a non- 
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equilibrium ^-entropy ifisl] . The successful application of Levy 
functions (TBWpp) to the spectra in p+p collisions at RHIC re- 
sulted in q values significantly larger than unity and are different 
between the groups of baryons and mesons 25]. However, 
in the central Au+Au collisions, the spectra at low pr show 
characteristic Boltzmann distribution with q value being close 
to unity — > 1, Eq.[T]becomes a Boltzmann distribution) even 
though there are still significant power-law tails with consider- 
able particle yields at high pr 1 10]. In addition to the escaping 
jets at high pr, coalescence with non-equilibrated quarks has 
also been proposed to study the power- law behavior . Dif- 
ficulty in accounting for these processes so far seems to be a 
major limitation of the TBW statistical description of the ex- 
perimental data over a wide pj range ifioil . To bridge the hy- 
drodynamic nature of the spectra at low pr and power-law tails 
at high pt with smooth transition at intermediate pr, models 
which include a hot and dense core with a corona of jet-like 
process have been proposed lfl4ll27ll . 



by extrapolating to unmeasured kinematic ranges. 
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Figure 1: (Color Online) Identified particle transverse mass spectra in p+p col- 
lisions (a) and 10-40% Au+Au collisions (b) at = 200 GeV. The symbols 
represent experimental data points. The curves represent the TBW fit. Only fits 
to the particles are shown since the model has the same spectral shapes for par- 
ticles and anti-particles. For plotting in panel (a), the spectra of meson (baryon) 
are scaled to match that of (p) at nix = 1.5 GeV/c~ for ^ (5-), at 4 GeV/c- 
for J /if/ and at 1 GeV/c~ for the rest. 

In this paper, we present the procedure of implementing 
nonextensive statistics in the Blast-Wave model (TBW) with 
azimuthally anisotropic particle emission, and use it to fit the 
identified particle spectra and for the first time to elliptic flow 
at mid-rapidity at RHIC. The model uses the TBW function ob- 
tained from p+p data ifioll as corona and an additional TBW 
function as core to fit Au-nAu data. The formalism thus pro- 
vides a systematic comparison between p+p and central A+A 
collisions in one macroscopic statistical model framework and 
gives an accurate numeric description of the experimental data 
over a wide range of pr- Examples of such successful applica- 
tions in the related subjects are the chemical fit to the particle 
yields 1 28 , 29[] and the global fit of the parton distribution func- 
tion (PDF) of proton ]30]. Good TBW fits can also offer a sim- 
ple formula for developing ideas and building models in a rea- 
sonably realistic environment flillllllllillllllllll, and 
provide a practical experimental tool to extract particle yields 
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Figure 2: (Color Online) Identified particle V2 in 10-40% Au+Au collisions. 
The X-axis is depicted by the kinetic energy (mr - m) 01, showing the scaUng 
of V2 at low Pt, and grouping of baryons and mesons at the intermediate pr 
range. The curves represent the TBW fit (Eq. [3j and the characteristic NCQ 
scaling. Also shown are the indistinguishable and A curves from TBWcore 
alone over the entire range. 

To take into account collective flow and azimuthal anisotropy 
in the transverse direction in relativistic heavy ion collisions. 
Levy distribution needs to be embedded in the framework of hy- 
drodynamic expansion ]20]. We follow the recipe of the Blast- 
Wave model provided in literature ]|5l S S H, and change 
sources of particle emission from a Boltzmann distribution to 
a Levy distribution ] 10]: 
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where yi, — \n(^Js^/m^) [35] is the beam rapidity and 
the rapidity distribution can be approximated as a Gaus- 
sian with a width of cr,, = 2.27 ± 0.02 at the center- 



of-mass energy of ^/s^ = 200 GeV |l36l |37|], transverse 
energy Et = mr cosh(3') cosh(p) - pr sinh(p) cos{(f)b - 4>), 
P = V('"cos(0s)//?x)^ + (rsin(0j)//?y)2(po + p2cos(20i,)) 
is the flow profile in transverse rapidity, and tan(0fc) = 
(7?x/^y)^ tan ((/I5) relates the azimuthal angle of the coordinate 
space {(f>s) to the angle of the flow direction {(pt,) of the emitting 
source Equati on [21 extends the nonextensive statistics in a 
blast-wave model flu] to incorporate particle emission from an 
elliptic source {Rx and Ry are the axes in the coordinate space) 
with an elliptic expansion (po and P2) |9|]. In addition to this 
core component, it is important to include the corona with jet- 
like particle emission at high momentum resembling an ensem- 
ble of individual p+p collisions The combined core 
and corona formula reads: 



dN 



mTdmjd(p ' 



= TBWr. 



+fppNhine{l + vf cos {2(l>))TBWpp, 



(3) 



where TBWcore is from Eq. |2] fpp and v^^' represent the frac- 
tion and the anisotropy of the escaping jet comparing to the 
expected number of binary p+p collisions in Au+Au collisions 
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Table 1: Values of parameters and best;^'^ from TBW fit of Eq.|?]to identified 
pailicle pj spectra in p+p collisions at RHIC. The uncorrelated systematic 
errors are included in the fit. 





T (MeV) 


q-l 




Mesons 


89.9 ± 0.7 


0.0955 ± 0.0006 


283/257 


Baryons 


68.5 + 4.0 


0.0855 ± 0.0019 


151/128 



(Nbin), e = PtUPt + 2^) takes into account the gluon satura- 
tion effect from p+p to Au+Au collisions with saturation scale 
2s = 1.5 GeV/c and TBW,,,, are the adopted fit results 
from the spectra in p+p collisions without any additional free 
parameter for either baryons or mesons. Since radial flow ve- 
locity p was found to be zero in p+p collisions jiotl . TBWpp 
can be simplified as: 



TBWpp cc niT \ dy e "V-^'t >'" cosh(y) 

J-y,, 

x(l + ^-^niT cosh(y))"'/<«"", 



(4) 



The STAR and PHENIX collaborations have published the 
most complete series of particle spectra and V2 at mid-rapidity 
for p-Hp and Au+Au collisions at -\fs^ - 200 GeV. The iden- 
tified particle spectra and V2 include tt'^, K"^, Ks, K*, p, (p, A, 
E,^ A, and H+ in STARpubHcations H B El, US El HI , 
ilSllillilliJIiilillsai. The T], Tj' and oj spectra in p+p 
collisions, K'^ spectra and 7r° V2 in Au+Au collisions are from 
PHENIX publications [SlllsllSEllilliS]. Figure[T]a shows 
the invariant differential yields in p+p coUisions at = 200 
GeV. All of the mesons or baryons have the common nij spectra 
shape. The solid lines represent the fit to Eq.|4]for mesons and 
baryons separately. The common fit parameters and hestx^ per 
fitting degree of freedom (nDoF) are listed in Tab. [T] In addi- 
tion to the common fit parameters, a global normalization factor 
was applied for each particle spectrum. The p^-integrated cross 
section dcr/dy for each particle species obtained from the fit are 
listed in Tab.|2] Figure[T]b shows the invariant differential yields 
together with our fit results in 10-40% centrality Au+Au colli- 
sions. The results of V2 from the simultaneous fit with spectra 
are displayed in Fig.|2] The fit parameters and best;^'^ per nDoF 
are tabulated in Tab. [3] In addition to these parameters common 
to all particles, a fit parameter is required as a normalization 
factor for TBWcore of each particle spectrum. The TBWpp in 
the fit function is directly adopted from the fit results shown in 
Fig.[T]a for each particle spices without any additional free pa- 
rameter. The Raa (ratio of the Nbi„ normalized pr spectra in 
A+A collisions to the underlying p+p spectrum) from model 
reproduces the data very well, as shown in Fig.|3] 

The Blast- Wave model with nonextensive statistics and az- 
imuthal anisotropy and a core-corona composition has allowed 
high quality fits to spectra and elliptic flow over a broad pr 
range. The striking feature of the experimental observations re- 
lated to the NCQ scaling is well reproduced. The results can 
be summarized as follows: 1) The bulk core alone (TBWcore) 
with finite q value can fit the data very well at low pr lflo[[l6ll . 
The new extension to V2 continues to provide high quality 
fits. The system produces maximum radial flow velocities of 



tanh(p() + P2) - 0.69c and tanh(po - Pi) - 0.43c along x-axis 
(the reaction plane) and y-axis, respectively. This demonstrates 
that the bulk system can be described with a few macroscopic 
parameters and is qualified as a thermodynamic state. 2) Fig- 
ure [3] shows that the Raa of the experimental data points at low 
Pt are below those at high pr- The e parameter necessarily 
brings the p+p component at low pr down to a subdominant 
fraction. This necessity may be partly due to the modification 
of the jet 10W-/97- component by the bulk fli]. 3) The non- 
equilibrated component in the bulk core produces a power-law 
tail in spectra and high V2 at the intermediate pr- 4) The baryon 
and meson yields {TBWpp) are grouped in p+p collisions with 
baryon yields systematically lower than meson yields as ob- 
served in the experimental data ifioi |27 , 41]. 5) The combina- 
tion of non-equilibrium tail from core and the baryon-meson 
separation from corona brings down the bulk V2, produces the 
baryon enhancement and the NCQ scaling at the intermediate 
Pt- 6) The medium quenches the jet and reduces it to a frac- 
tion (Jpp - 0.36) of its underlying binary nucleon-nucleon col- 
lisions, resulting in a finite azimuthally anisotropic emission 



{v'°' = 8.7%) 1601]. 7) The high-precision (±10-"* stat.) ex- 
perimental data points concentrate at low pj, dominating the 
fit x^- Additional high-quality data in higher pj range (e.g. 
V2 of baryons) will balance the contributions to the x^ from 
components with different physical origins. The not-quite-ideal 
X^lviT)o¥ value indicates significant tensions among the differ- 
ent datasets and model, and warrants further detailed assess- 
ment and categorization similar to the PDF fit in the future. 

This TBW model describes the number of constituent quark 
(NCQ) dependencies observed in the data which have also been 
described in recombination models lH , 12 , 2^ 62 ] . The NCQ 
dependence indicates that baryon production at intermediate Pt 
increases faster going from p+p to Au+Au collisions than me- 
son production does. The division into a core and corona com- 
ponent in Eq. |3]represents a division into two samples: a core 
and corona where the core has a larger fraction of baryons than 
the corona. We describe both the core and corona with a TBW 
statistical model. But whereas the recombination picture at- 
tempts to provide a microscopic description for why the bulk in 
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Figure 3: (Color Online) Comparison of TBW fit (lines) with STAR measure- 
ments of Raa in 10-40% Au+Au collisions. The syst. errors (horizontal bars) 
and Stat, errors (vertical lines) are shown separately for protons and charged 
pions. Systematic and statistical errors are added quadratically in the fit, and 
are shown as such in the plot except the protons and charged pions. 
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Table 2: Cross sections dtrldy in mb of different particles at mid-rapidity in p+p collisions at yfs^ = 200 GeV. "TBWpp" column show the values obtained 
from fit results shown in Fig.^a. "Data Ref." column shows where are the data points from. "PHENIX" column shows the values obtained by fits to PHENIX 
measured spectra with the Tsallis functional form as described in Ref. |56]. "Published" column shows the values published by experiments. An additional 9.7% 
systematic uncertainty should be added to all dcr/dy values listed in the column "PHENIX" and 12% to the values in the column "'STAR' to account for the 
trigger uncertainties. Values in the column "Pubhshed" are also given without these systematic uncertainties. The column "S.M." is the prediction of the statistical 
model fSTt]. All errors are the combined statistical and systematic uncertainties. 



Particle 




Data Ref. 


PHENIX [56] 


STAR 


S.M. [57] 


tt" 






41.4 + 5.8 




46.9 


n"- 


42.4 


[5,47] 


39.4 + 7.3 


43.2 + 3.3 [34] 


42.1 


7r~ 


41.9 




38.6 + 7.2 


42.6 ±3.3 [34] 


41.5 


K+ 


4.69 


[5,44,413 


4.57 + 0.61 


4.50 + 0.39 [34] 


4.57 


K- 


4.57 


[5,44,413 


4.20 + 0.51 


4.35 + 0.39 [34] 


4.38 


K° 

s 


4.73 


[411 


5.28 + 0.53 


4.02 + 0.34 [41] 


4.40 


V 


4.94 


[51] 


4.47 + 0.96 




4.93 


p 


6.79 


[501 




7.8+ \.2[50} 


5.58 


0) 


3.75 


[52j 


3.65 + 0.77 




5.03 




0.58 


[56] 


0.62 + 0.17 




0.365 




1.51 


[39] 






1.57 




1.54 


[39] 




1.52 + 0.19 [39] 


1.55 




0.521 


[45] 


0.421 +0.055 


0.540 + 0.086 [45] 


0.339 


7/(A(xlO^) 


0.741 


[53] 


0.759 + 0.053 






i//' (xlO^) 






0.133 + 0.031 






P 


3.96 


[5,47] 




4.14 + 0.30 [34] 


4.47 


P 


3.19 


[5,47] 




3.39 + 0.36 [34] 


3.59 


A 


1.36 


[41J 




1.31+0.12 [41] 


1.30 


A 


1.21 


[41J 




1.19 + 0.11 [41] 


1.11 




0.098 


[41J 




0.078 + 0.028 [41] 


0.092 




0.094 


[41J 




0.087 + 0.031 [41] 


0.082 


2*+ + 2*- 


0.359 


\M 




0.321 +0.044 [58] 


0.308 


£•+ + f- 


0.313 






0.267 + 0.038 [58] 


0.260 


A* + A* 


0.125 






0.104 + 0.017 [58] 


0.168 


Q- (xlO^) 


12.3 


[41] 




10.2 + 5.7 [41] 


17.1 



Table 3: Values of parameters and best;^'^ from TBW fit of Eq.[3]to identified particle pr spectra and V2 in 10-40% Au+Au collisions at RHIC. The uncorrected 
systematic errors are included in the fit. Syst. error of STAR V2|£P) are taken to be 15%/ Vl2 148]. The spectra contribute 359/244 to the;f2/nDoR Qs = 2.1 ± 0.2 
GeV/c with;>^^/nDoF=497/295 when Qs is a free parameter Including the V2IEP] of K- and p at low pr 159] increases the best;t^/nDoF to 769/307. 



PO 


P2 


Rx/Ry 


T (MeV) 


q - 1 fpp 


v^' ;^'^/nDoF 


0.654 + 0.002 


0.199 + 0.002 


0.871 +0.004 


128 + 2 


0.044 + 0.001 0.36 + 0.01 


(8.7 + 0.4)% 506/296 



Au+Au collisions has a larger fraction of baryons, the macro- 
scopic TBW statistical model only describes the existing states 
as they are, and is not sensitive to the underlying mechanism re- 
sponsible for producing those states. Regardless of whether the 
TBWpp component is indeed from individual p+p collision or 
part of the hadrons from coalescence of the escaping partons, 
Eq. [3] does provide a unified tool to describe simultaneously 
a variety of measurements over a wide range in pr for future 
studies. These results also provide a macroscopic foundation 
for discussing the entropy issues associated with the system's 
underlying microscopic subprocesses 123, 2^ 631,641 . 

In summary, we have implemented the nonextensive statis- 
tics in a Blast-Wave model and incorporated a core-corona 
model to describe a complete data set on identified particle 
spectra and elliptic flow versus transverse momenta at mid- 
rapidity measured at RHIC. Such a formalism simultaneously 
describes several novel observations reported at RHIC over a 



broad pr range. Specifically it provides an alternative phys- 
ical picture, through the role of the core-corona and non- 
equilibrium effects in understanding the baryon-meson differ- 
ences in the nuclear modification factor and the elliptic flow 
at intermediate pr- The macroscopic nonextensive statistics 
provides a complementary microscopic method to the hydro- 
dynamics and hadronic cascade for studying the evolution of 
the nucleus-nucleus collisions and the properties of QGP 
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